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Abstract. We consider the focusing nonlinear Schrodinger equation, 
in the L^-critical and supercritical cases. We investigate numerically 
the dependence of the blow-up time on a parameter in three cases: de- 
pendence upon the coupling constant, when the initial data are fixed; 
dependence upon the strength of a quadratic oscillation in the initial 
data when the equation and the initial profile are fixed; finally, depen- 
dence upon a damping factor when the initial data are fixed. It turns 
out that in most situations monotonicity in the evolution of the blow-up 
time does not occur. In the case of quadratic oscillations in the ini- 
tial data, with critical nonlinearity, monotonicity holds; this is proven 
analytically. 



1. Introduction 

Consider the Schrodinger equation with focusing nonhnearity {X,a > 0): 

(1.1) idtu + Au = -X\u\'^''u , {t,x) £R+ xW^ ; uit=o = uq . 

Such an equation may appear as an envelope equation in the propagation of 
lasers (see e.g. and ^lEl for a rigorous mathematical justification). 
It is well known that if uq G i7^(R") (inhomogeneous Sobolev space) and 
a < then ()1.1() has a unique solution in H^{W^), defined locally in 

time (see e.g. jSj)- It needs not remain in H^{W^) globally in time: finite 
time blow-up may occur when a > ^ and A > (focusing, or attractive, 
nonlinearity); see e.g. El EHI- Since A S M, the L^-norm of n(t, •) is 
independent of time, and finite time blow-up means that there exists T > 
(we consider only forward time evolution) such that: 

ll^a:'u(t)||L2 — > -|-00 aS t ^ T . 

Many papers provide important properties about the blow-up rate or the 
profile at blow-up; see e.g. HHSOlIIHlIIllIZlllllilinilllllinifor some an- 
alytical results, and fo'^ ^ i^i^e survey of the latest results. We also point 
out that recent numerical experiments T21 have shown a completely new 
phenomenon, where the blow-up profile is independent of the usual ground 
state (quasi self-similar ring profile), and the blow-up rate seems to be the 
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minimal one given by the theory (square root blow-up rate). As recalled in 
j27j . the main three directions of research in this subject are: giving suffi- 
cient conditions to have finite time blow-up in the energy space; estimating 
the blow-up rate and the stability of the blow-up regimes; describing the 
spatial structure of the singularity formation. On the other hand, it seems 
little attention has been paid to the time where blow-up occurs. In this pa- 
per, we investigate by numerical experiments the dependence of the blow-up 
time upon, for instance, the coupling constant A, when the initial datum mq 
is fixed. 

To motivate our study, we recall some results from and . In [H] , the 
authors prove that if the initial datum uq{x) is replaced by tio(x)e~*^l^l 
then one can relate explicitly the blow-up time of the corresponding new 
solution Uf, to that of u, in the case of a critical nonlinearity, a = ^. This is 
a consequence of the conformal invariance. In the super-critical case a > 
the conformal transform does not leave invariant, and introduces a 

factor (1 — ht)"''^~'^ in front of the nonlinearity. It is also established in |5] 
that if u has negative energy (in this case, there is finite time blow-up at 
least if xuq G L^(R"') jTl]; see also 123 , 17^ for weaker assumptions), then 
for large b, blow-up occurs sooner than for 6 = 0; unlike in the conformally 
invariant case, one does not know whether the blow-up time is monotonous 
with respect to b. The numerical experiments we present here show that it 
is not monotonous with respect to b. 

In the author considers the damped cubic Schrodinger in space di- 
mension two: 

(1.2) idti^ + Ai; = -id^j - 1-01 V , {t, x) S M+ x ; ip\t=Q = uq . 

It is conjectured that the blow-up time is monotonous with respect to (5 > 0, 
and guessed that the same holds in the super-critical case. Our numerical 
experiments show that neither of the two guesses is satisfied. These two 
guesses are very satisfactory when one think of the initial datum uq as a 
single hump; numerics also suggest that in this case, the blow-up time is 
monotonous with respect to 5. On the other hand if uq is made of, say, two 
humps, then the initial intuition seems to be wrong. 

Note that introducing u{t, x) = e^^'ip{t, x), Equation (|1.2j) is equivalent to: 

(1.3) idtu + Au = -e~^^^\u\'^u , {t,x) £R+ xR^ ; uit=o = uq . 

As in the case of initial quadratic oscillations mentioned above, this trans- 
form yields an equation of the form (|1.1|) . with a time-dependent coupling 
"constant", A = e~^^*. This function of time is monotonous, decreasing. 

This naturally leads us to the question, when A is really a constant: is the 
blow-up time monotonous with respect to A? Tracking the dependence of 
the blow-up time upon A can be viewed as both a generalization and a study 
case of the two problems raised in [5] and Numerics show that both in 
the critical (c = ^) and in the supercritical case {a > one should not 
expect the blow-up time to be monotonous with respect to A. Essentially, 
the idea is the same as what was announced above: when the initial datum 
Uq is a single Gaussian, then monotonicity seems to holds. When uq is the 
superposition of two such functions, there is a lack of monotonicity. 
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Roughly speaking, suppose that two humps are placed to both sides of 
the origin and that the natural (free) time evolution tends to send mass 
from both these humps to the origin. This may occur thanks to a particular 
phase term, for instance, or simply mass dispersion. For large A, the blow-up 
occurs before the two humps have merged: the two humps do not interact 
before the blow-up, and break down independently. For small A, these humps 
merge into one hump before blow-up takes place. Asymptotically, for "large" 
A or for "small" A, we observe some monotonicity in the blow-up time. On 
the other hand, the monotonicity breaks down in the "transition" region, 
that is, for intermediate A. 

The rest of this paper is organized as follows. In Section [2 we recall 
some analytical results on local existence and finite time blow-up for the 
nonlinear Schrodinger equations that appear in this paper. For the sake of 
readability, proofs are given in an appendix. Numerical tests are presented 
in Section |21 to [3 In Section |31 we consider the dependence of the blow-up 
time upon A in in Section^ we measure the dependence of the blow- 

up time upon a modification of the initial datum with quadratic oscillations; 
Section is devoted to tests on l\l.'A^ and its natural generalization in one 
space dimension. 

2. Some theoretical results 

In this section, we present some analytical results that provide bounds, 
from above and/or from below, for blow-up time. The techniques are clas- 
sical: the proofs rely on Strichartz estimates, conservations of mass and en- 
ergy, and pseudo-conformal conservation law. Yet, it seems that the explicit 
dependence of the existence time upon some parameters had not been inves- 
tigated before, except in [Hj (see Section . These results are somehow a 
more quantitative motivation for the numerical tests that follow. Technical 
proofs are given in Appendix 1X1 

2.1. Generalities. We first state a local existence result, from which we 
infer a lower bound on the blow-up time. Recall that in we consider 

only nonnegative time. 

Proposition 2.1. Let X > 0, a > ^ with a < if n > 3, and uq G 
H^{W^). There exists e > independent of X,a and uq such that if 

2-(n-2)a _i^(2-(n-2)cr) ,, „^ 

(2.1) XT^^WuoWlt'' llV^noll^r <£, 

then (HHJ) has a unique solution u G C {[f),T[]H^) f] {[0,T[; 1'^'^+^) . 

If moreover 

noGS:={/GFi(M") ; \x\f e L\R"')} , 

then u £ C ([0,T[; S). In addition, the following quantities are independent 
of time: 

(2.2) Mass: M = ||n(t)||i2 = Const = ||no||L2, 

(2.3) Energy: E = \\V Mt)\\l2 " -^^\\u{t)\\%tl2 = Const. 
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This result is standard (see e.g. "S"), except that usuahy Condition (12.1(1 is 
not given exphcitly: see Section FA.ll for the proof. Note that in Corollarv l2.2l 
below, we modify this condition, to a somehow weaker one. The reason is 
that in Proposition 12.11 local solutions are constructed without taking the 
conservations of mass and energy into account. We shall not recall explicitly 
the fact that for small initial data, the solution to p.l|l does not blow up. 
Indeed, for small data, the conservations of mass and energy yield an a priori 
bound on the H^-norm of the solution, thus ruling out finite time blow-up. 
In our case, small initial data means that, for instance, A-'^/^^'^^ ||uo||/fi < 6 
for some constant 6 depending only on n and a. We shall not insist on that 
aspect, since our analysis is focused on regimes where blow-up does occur. 
This is why in the next two corollaries, A is morally "large", while uq is 
fixed. 

Corollary 2.2. Let X> 0, a > ^ with a < if n > 3, and uq £ H^{W). 
There exists e > independent of A such that if 

(2.4) (A^'^llV^nolli'^ + 1) r2-("-2)- < e , 

then (dl} has a unique solution u G C {[0,T[;H^) n {[0,T[; 1'^"+'^) . 

If moreover uq E S, then u £ C ([0, T[; S). 

The proof of Corollary 12.21 is given in Section rOl 

Corollary 2.3 (Dependence with respect to the coupling constant). Let 
A > 0, (7 > ^ with a < ifn>3, and uq G S. Assume that u blows up 
in finite time T* > 0. 

2n- 

1. We have T* > C (A) 2-(n-2)<T^ jgj- gome constant C independent of X, 
where (A) = \/l + A^. 

2. If in addition E < 0, then T* < C (A)"^/^ for some constant C inde- 
pendent of A. 

Remark. The above constants C and C are independent of A, but depend 
on the other parameters, uq, n and a. As a matter of fact, the construction 
we use yields constants depending only on ||uo||s! n- and a. 

Proof. The first part is a straightforward consequence of Corollarv l2.2l The 
second follows from the Zakharov-Glassey method |32| I14j . Introduce 

y{t) = / |a;p|u(t, x)|^(ix . 

Then since A > and o" > -, we have: y{t) < 8naE, where E denotes 
the energy defined in Proposition [23 Integrating, we infer that for large A, 
y(t) < —CXt^, for some positive C independent of A. Since y{t) > so long 
as u remains in S, this yields the second point of the corollary. □ 

We always have 2-{n^2)a > ^) so the above two bounds go to zero with 
different rates when A +oo. Condition (|2.1() would yield only T* > 

__2ct+2_ 

C (A) '2-{n-2)a _ "Wlthout cvcu trying to see if any of these bounds is sharp, 
we ask the following question: 

Question 1. For a > - and a fixed initial datum uq G is the blow-up 
time for u solution to ()1.1() monotonous with respect to X? 
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This issue is addressed numericahy in Section |31 where our resuhs show 
that the answer to the above question should be no. 

2.2. Initial data with quadratic oscillations. Like in P|, we now fix the 

equation, and alter only the initial data, with quadratic oscillations: 

(2.5) idtu + Au = -\u\'^'^u , X €R"- ; uit=o = uoix) , 
with no G S. For a / 0, define: 

/ \ / \ e^^t-a) / at X \ a — t 

(2.6) vlt, xj = —— — j-u 'TTT ) where h(t) = • 

^ ^ ^ ^ /i(t)"/2 \a-t h{t) J ^ ' a 

Then v solves: 

(2.7) idtv + Av = -h{t)''''-^\v\'^''v ; vit=o = uo{x)e-'^ . 

In the conformally invariant case a = 2/n, v solves the same equation as u. 
The only difference is the presence of (additional) quadratic oscillations in 
the data. 

Proposition 2.4. Let € S and 2/n < cr < 2/(n — 2). Suppose that u 
blows up at time T > 0. Let a G M*. 

• If a> 0, then v blows up at Ta(v) = —T < T. 

a + T 

• If a < and a + T < 0, then v blows up at Ta(v) = — T > T. 

a + T 

• If a < and a + T > 0, then v is globally defined in S for positive 
times (but blows up in the past if a + T > 0). 

For the critical case a = 2/n, this result is proved in [Hj (see also 0). We 
sketch a slightly different proof in Appendix I A. 31 which easily includes the 
case 2/n < a < 2/(n - 2). 

In the super-critical natural question is to understand the role of 

the function h. Introduce w solving: 

9 _ -Mi 

(2.8) idtw + Aw = -\w\ w ; w\f=Q = uo{x)e * . 

Proposition 2.5. Let G S and 2/n < a < 2/(n — 2). 

• Let Ta{w) denote the maximal existence time in the future for w. 
Then there exists C = C{\\uQ\\Y,,n, a) independent of a such that: 

Ta{w) > C|a|2~("-2)-. 

• // the energy E of u is negative, then for a > 0, w blows up at time 
Ta{w) < a. 

The first point is a direct consequence of Corollarv l2.2l The second point 
is proven in [9,,, and relies on the pseudo-conformal law for w. 

To understand the infiuence of the quadratic oscillations on the blow-up 
time, we have to compare the blow-up time of u and that of w. In the critical 
case, the blow-up time depends explicitly on the magnitude of the quadratic 
oscillations via Proposition 12. 4( since f = it; by conformal invariance. In the 
super-critical case, we ask: 

Question 2. For c > f and a fixed uq G S, is the blow-up time for w 
solving 1)2. 8|) monotonous with respect to a? 
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This issue is addressed numerically in Section |11 we first compare the 
numerics with the analytical results in the conformally invariant case, then 
perform tests in the supercritical case which indicate that the answer to the 
question above should be no. 

2.3. Damped equation. We now consider: 

(2.9) idfil^ + = -l-ipf^tp -iSip , (t,x)GM+xM" ; '0|t=o = ^^o , 

for 5 > 0. A direct application of (the proof of) Proposition 12. II shows that 
tp is defined on [0,r[ for T > C\5\-^. Note that the sign of 6 is irrelevant at 
this stage. To take damping effects into account, introduce 

u{t,x) = e^^Tp{t,x) . 

Then u solves: 

(2.10) + Au = -e"2<T5t . uit^Q = uo. 

For 5 sufficiently large, u is defined globally in time, in the future: 

Proposition 2.6. Let <7 > | with a < if n > 3, and uq € H'^{R"-). 
There exists C = C{a,n) depending only on a and n such that if 

(2.11) 5"+' > C\\uof^t^''-^^'''^\\Vuo\\lf, 

then fnn^ has a unique solution u€ C (M+; H^) Pi (M+; 1'^"+'^) . 

The proof is given in Section I A. 41 The following question is addressed 
numerically in Section |S1 

Question 3. For a > - and a fixed uq G S, is the blow-up time for u 
solution to (|2.1U|) monotonous with respect to 6 > 

The simulations show that the answer is no. 

3. Numerical Test, dependence on A 

We perform numerical tests by using a direct discretization method for 
equation (jLlj) . respectively (|2.1U|) . By this approach no restrictions are im- 
posed on simulations on the closeness to eventual blow-up points in space 
and time and by the evolution in places "not close" to the blow-up. We em- 
ploy two different numerical methods: the Time-Splitting Spectral method 
(TSSP), and the Relaxation method (RS). 

The TSSP is based on an operator splitting method, the split-step method. 
The split-step method is based on a decomposition of the flow of the non- 
linear equation (|l.lj) (or (|2.1fljl ). Define the flow as the flow of the linear 
Schrodinger equation 

f idtv + Au = 0, xeW, t>Q, 
\ v{x,Q) = vo{x), X G M^, 

and y* as the flow of the nonlinear differential equation 

r idtw = -g{t)\w\^''w, xGM", t>0, 
\ w{x,Q) = wq{x), X e M^. 

where g{t) = A for the case of (mU) and g{t) = e'"^"^^ for (tTTTH) . Then the 
split-step method consists of approximating the exact u{x,t) at each time 
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step by combining the two flows X* and Y^. We employ here the Strang 
formula = which is of second order. Higher order splitting is 

also possible. For a general nonlinear Schrodinger equation, a convergence 
proof of the splitting method was done by Besse et al. in For the 
TSSP, a spectral method is employed to compute the flow of the free 
Schrodinger equation. The flow Y^, which is the flow of a nonlinear ODE, 
can be computed exactly, since it leaves \w\'^ invariant. So its integration is 
straightforward. The TSSP has proved to be an efficient and reliable method 
for NLS type equations. See for example for a study of the NLS in the 
semi-classical limit case, and j251 [Sj for a more general numerical study. 

For the 2-d calculations, a parallel version of the TSSP scheme is used on 
the parallel cluster machine "Schrodinger III" at the University of Vienna. 

The Relaxation method (RS) is a discretization of finite difference type 
[HI. It is based on central-difference approximation shifted by a half time- 
step. To formulate the scheme, we rewrite the NLS equation as a system: 

idtu + Au = —g{t) V' u, 

where g{t) is as above. Then the first equation is discretized by central 
difference quotient at the time tn+i/2 = {n + ^)At, and the second by 
central time average at the time tn = nAt. Let be the approximation at 
t = tn, then the scheme is given by 

2 ' 

where Ad denotes a finite difference Laplacian. It requires only explicit 
evaluations of the nonlinear term. It also conserves energy ( ) . 

To determine whether blow-up is occurring or not, we calculate the two 
terms in the energy (|2.3)) . kinetic and potential energy, and look for an 
increase of at least four orders of magnitude in both of them. The first time 
at which this is occurring is assumed to be the blow-up time. The space and 
time resolution are chosen sufficiently fine such that this increase is realized. 

3.1. Critical power. First we consider (jl.lj) for cr = ^, that is the critical 
case. We study the dependence of the blow-up time on the constant A for a 
series of different data. 

Tests in one space dimension 

Test 1. For the case n = 1, the first kind of data we study is 
(3.1) no(x) = Ce"^' g-ilog(e-+e-)^ 

The constant C is equal to 1.75 which leads to ||Mo|li2 = 3.839. We take 
Np = 2^^ = 4096 mesh points, and several time steps according to the 
necessary resolution for the blow-up, with At = 2.5 • 10~^ as smallest. The 
discretization domain is [—8, 8]. Figure^shows the blow-up time in relation 
to a changing A. It can be observed that the blow-up time is decreasing 
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critical power 

0-16 I 1 1 1 1 




Figure 1. Blow-up time with varying A, single Gaussian 
data (Test 1). 



critical NLS, two hump data 

0.45 1 1 1 1 1 1 

0.4 - 




0.1 I ' ' ' ' ' ' ' ' 1 

0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 

Figure 2. Blow-up time with varying A (Test 2). 

monotonously with A, as predicted by the heuristics for the case of a single 
Gaussian profile. 

Test 2. The next kind of data we study is 
(3.2) uo{x) = C (e-^' - 0.9e-3^') g-iiog(e-+e--)_ 

The constant C is equal to 4 which leads to ||tio||^2 = 3.907. The difference 
of two Gaussian profiles results in two local maxima in the modulus of uq. 
The phase term has a focusing effect and its focus point does not agree with 
the local maxima of the modulus. The finest discretization parameters used 
are Np = 2^^ = 16384, At = 2.5 • IQ-^. The discretization domain is [-8, 8]. 

Figure [21 shows the blow-up time in relation to a changing A. It can be 
observed that for low and very high strengths of the nonlinearity A (A < 1.6 
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Figure 3. Time evolution with low potential energy. 



S, . . . _ ; " ' ' " .. L ' ■ • 




Figure 4. Time evolution with high potential energy. 

and A > 2.2), the blow-up time is monotonously decreasing with A, while in 
between there is a region where monotonicity does not hold. 

Heuristically, two effects play a role in the solution with this data: the 
nonlinear self-focusing, which tries to focus the mass to points where the 
most mass is already present, and the (linear) phase influence, which tends 
to focus the mass at zero. 

In Figures 13 Sand [51 the time evolution of the modulus of u{t, x) is shown 
for values of A from the three different regions of the above curve. 

• Figure 131 shows the case A = 1. The two initial humps merge to one hump 
before the blow-up, which happens at a single point. The phase focusing 
happens at a faster time scale than the nonlinear focusing. 

• Figure 0] shows the case A = 2.4. Blow-up is occurring simultaneously 
at two points and there are always two humps present. Apparently the 




Figure 5. Time evolution in intermediate regime: non-monotonicity. 



X = 2.0 

140 



120 



100 
80 
60 
40 
20 



'0 0.05 0.1 0.15 0.2 0.25 

t 

Figure 6. Determining blow-up time by increase of energy: 
(minus) total energy in comparison with kinetic and potential 
energy away from the blow-up 

nonlinear self-focusing here happens faster than phase focusing. 
• Figure ini is for A = 2.0, which is in the non-monotonicity region. Blow-up 
here occurs at a single point, and the merging of the two humps is closer to 
the blow-up than in Figure El In this case it is not clear which of the two 
effects would happen at a faster time scale, nor how they would interact. 
Blow-up is occurring, but the blow-up time is no longer monotonous with 
respect to the size of the nonlinear term. 

In Figure ini we show the energy components for the case A = 2.0 which 
are used to determine the blow-up time. The kinetic energy increases by a 



kinetic energy 
potential energy 
- (neg.) total energy 
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critical NLS, ttiree iiump data 



0.4 



0.35 - 



0.2 - 



0.05 



0.15 - 



0.25 - 



0.3 - 
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0.5 



.5 



2.5 



3.5 



Figure 7. Blow-up time with varying A, three hump data (Test 3). 

factor of 10'^, or shghtly more, at the blow-up time. Figure El compares the 
kinetic and potential energy parts with —E[u{t)) away from the blow-up. 
The TSSP scheme used for this simulation does not conserve energy. Near 
the blow-up time, the energy conservation is no longer true for the numerical 
result. Note that, outside of a time interval close to blow-up, the kinetic 
energy is not monotonously increasing in time. 

Remark. If the data (|3.2j) are used without the phase term, which leaves just 
a "two-hump" profile, non-monotonicity can be observed in the same way as 
described above. However overall blow-up times increase. Apparently the 
merging of the two humps can occur even without the influence of an initial 
phase term thanks to the mass dispersion tendency of the free evolution 
together with the focusing effect of the nonlinearity. 

Test 3. Data with three humps 
The next test uses a sum of three Gaussians and the same phase term as 
before, so there are three local modulus maxima instead of two. 



with C = 2, ||uo||^2 = 5.09. The discretization parameters are Np = 2^"^ 
mesh points and At = 1.5 • 10~^. The discretization domain is [—8,8]. The 
blow-up times with respect to changing A are shown in Figure [3 The same 
effect as above can be observed. 

There are two regimes for the nonlinearity strength A where the blow-up 
time is non-monotonous. Blow-up happens either at three points, two or 
one point depending on A. 

Test 4. Data with two humps up, one down 
In this test the data are taken to be 
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critical NLS, ttiree iiumps two up, one down 
0.35 1 , , , , r 



0.3- 




1.6 1.8 2 2.2 2.4 2.6 2.8 3 

X 

Figure 8. Blow-up time with varying A, three hump data (Test 4). 

with C = 2, ||uo||^2 = 4.93. One of the three Gaussians has an opposite 
sign, so there is a constant phase shift in part of the data. The blow-up 
times are shown in Figure |H1 

For A < 1.6, there is no blow-up occurring. For larger A, non-monotonicity 
similar to the situation above can be observed. Observe that the slope of 
the curve is rather steep, for 2.7 < A < 2.8. This shows a highly nonlinear 
phenomenon. Note however that there is one point on the steep line after 
the local maximum: 



A 




2.7 


0.1912 


2.725 


0.2152 


2.75 


0.1612 


2.8 


0.0555 



Test 5. Data with one hump up, one down 
We use 

uo{x) = Ce"^' tanhxe-^'°s(e'+""'). 
We use C = 3.0 and ||uo||^2 = 4.4476, and discretizations of Np = 2^^ mesh 
points and At = 2.0 • 10~^. The result is shown in Figure 1^1 In this case, 
the blow-up time is monotonous with A. 

Remark. For other point-symmetric data with more than two humps, there 
is monotonicity. 

Test 6. Stability issues: asymmetric data 
In order to investigate the dependence of the qualitative findings on the 
symmetry of the data, we choose some data which perturb the symmetry of 
(lOl . We used 

Uo{x) = C (e"^' - 0.9e-3^') g-ilog2cosh(x-0.25)^ 

and 

no(x) = C (e-(^-i-5)' +0.99e-(^+^-5)') ^-.iog(e^+e-)^ 



BLOW-UP TIME FOR NLS 



13 
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Figure 9. Blow-up time with varying A, data with point 
symmetry (Test 5). 

to have either a phase that focuses away from the symmetry point between 
the two humps, or two humps of different heights. It turns out that non- 
monotonicity can be observed for for both of these cases, too. The pertur- 
bations to (|3.2j) tested here are smah enough not to hinder non-monotonous 
blow-up times. 

Test in two space dimensions 

Test 7. For the test in two space dimensions, we extend (|3.2j) by making 
the phase term radially symmetric and multiplying the one-dimensional two- 
hump profile by a single Gaussian in the second space dimension: 

(3.3) uo{x,y) = C (e--' - 0.9e-^^') ^-,2^-.iog2cosh(V^) ^ 

Figure [TUl shows the modulus of H3.3|) . We choose C = 7.0, ||n|||2 = 15.0. 
The smallest discretization parameters used are Np = 2^^ mesh points and 
At = 1 • 10~^ with the discretization domain [—4,4]^. 

The blow-up times with changing A are shown in Figure lTTl Non-monotonicity 
can be observed. 

3.2. Supercritical power. 

Test 8. We tested equation in one space dimension, with cr = 3 

and the data (|3.2j) {C = 3.5, hence ||tto||^2 = 2.99). The discretization 
parameters are Ax = 0.0039, and up to At = 1 • 10~^. The discretization 
domain is [—8,8]. The blow-up time with varying A is shown in Figure IT21 
Also in the supercritical case, non-monotonicity of blow-up times can be 
observed. 

4. Numerical Test, dependence on quadratic oscillations 

We turn to equation 1)2 .71) and investigate the dependence of blow-up 
time on the scale of quadratic oscillations. To compare the simulations to 
the result of Proposition 1^31 we simulate H2.5() with the same data to obtain 
the blow-up time for this equation, and then plot the curve of Tq that is 
predicted in Proposition 12.41 for the two regions a > and a < 0, a + T < 0. 
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Figure 10. Initial data for two dimensional case, modulus (Test 7). 



critical two-dim. NLS 
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Figure 11. Blow-up time with varying A, two dimensional 
case (Test 7). 



4.1. Critical power. We use (12. 7p in space dimension one, with a = 2 
with various uo{x). 

Test 9. Single hump: We take 



uo{x) = C e 
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supercritical NLS 

0.35 I , , , , , 




4 0.6 0.8 1 1.2 1 4 1.6 1.8 2 2.2 2.4 

Figure 12. Blow-up time with varying constant A (Test 8). 

with C = 1.75. The discretization parameters in this and the two following 
tests are Np = 2^^ mesh points and At = 4 • 10~^. The discretization 
domain is [—8, 8] for all cases except the two largest negative a in Test 9 
and 10. Here the domain is extended to [—18, 18] and the space resolution 
is Np = 2^^ resp. [—40, 40] and Np = 2^^ for the largest negative o of 
Test 11. Figure IT^ shows the blow-up time of v in relation to the scale a 
of quadratic oscillations in the data. We use both positive a and negative 
a with a + T < 0. Asterisks denote the blow-up times and the dashed line 
shows the result of Proposition 12. 41 with T obtained by a simulation of ()2.5|) . 
It can be observed that the results agree very well. In each of the two regions 
for a that have been used in this test, the blow-up time is monotonous with 
respect to a. 

Test 10. Two humps: We take 

uo{x) = C (e"^" -0.9e~3x2^ 

which is the same as 1)3. 2() without the phase term that appears there. We 
take C = 4.0. Figure shows the blow-up times marked by asterisks. 
As above the dashed line shows the modified blow-up time according to 
Proposition 12.41 The two curves agree. 

Test 11. Two humps with additional phase: Here we use 

uo{x) = C (e-^' - 0.9e-3^') g-iiog(e-+e--) 

which is the same as (|;-{.2j) . Figure El shows the blow-up times of the simu- 
lations and according to the Proposition as above. The results agree. 

4.2. Supercritical power. We now consider Eqns. 1)2. 7|) and (|2.8() with 
(7 = 3, in space dimension one. We give results for both w and v. We use 
the same series of data as for the critical power case. 

The discretization parameters in Tests 12 to 14 are Np = 2^^ mesh points 
and At = l- 10'^, down to At = 4 • 10"^. 

Test 12. Single hump: We take uq{x) = C . Figure [TBI shows the 
blow-up time in relation to the scale a of quadratic oscillations in the data. 
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0.45 
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Figure 13. Blow-up time with varying a in quadratic os- 
cillations, critical power (Test 9). 



3.5 



3 - 




Figure 14. Blow-up time with varying a in quadratic os- 
cillations, critical power, two hump data (Test 10). 



In the left figure, the dashed line is the calculated blow-up time for v and 
the asterisks mark the simulated blow-up times. The right figure compares 
the blow-up times for v and w, where w is denoted by asterisks and v by 
dots joined by a line. We can see that w blows up a bit earlier than v for 
positive a. For negative a, it blows up a bit later, but the blow-up times are 
rather close to those of v as in the positive a case. Again the blow-up times 
for V match the result of Proposition 12.41 

Test 13. Two humps: We take ^0(2;) = C (^e"^" - 0.9e-3^') . Figure HH 

shows the blow-up time of v in relation to a. 

Test 14. Two humps with additional phase: Here we use 

uo{x) = C (e-"' - 0.9e-3-') e-*i°g(-^+^-^). 

Figure ITHl shows the blow-up time of v in relation to a. 
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Figure 15. Blow-up time with varying a in quadratic os- 
cillations, critical power, data with additional phase 
(Test 11). 

Test 15. Two humps placed asymmetrically: Here we use 
«o(:c) = c(e-(3-)^+e-(3(-i-5))^]. 

We take C = 1.8. The smallest discretization parameters used in this test 
are Np = 2^^ mesh points and At = 1.5 • 10~^. The computation domain is 
[—8,8], except for the largest value of a. Figure IT^ shows the blow-up time 
for w in relation to a obtained by the two numerical methods employed: the 
circles represent simulations done by the TSSP, the asterisks simulations 
by the RS. In addition a solid line displays the blow-up times for v. Non- 
monotonicity can be observed, which answers question 2 in a negative way. 
The blow-up time of w is always smaller than that of v. We also see that 
the results of the two different schemes agree in a good way. Note that the 
occurrence of non-monotonicity is very sensitive to the size of the data. If 
we choose C = 1.7 or 1.9 instead of C = 1.8, monotonicity can be observed. 

Remark 4.1 (Test 15). Leave out the question of quadratic oscillations, and 
consider u solving 1)2. 5() . Its blow-up time as a function of the constant C is 
given by: 



c 




1.795 


0.528 


1.798 


0.480 


1.8 


0.462 


1.804 


0.446 


1.808 


0.507 


1.81 


0.076 


1.82 


0.048 



Since changing C in 1)2. 5() is equivalent to changing A in ()1.1() . we see here 
a behavior analogous to Test 8 (Figure IT^ . where similar data are used. 
The value C = 1.8 is very close to (actually slightly below) the potential 



18 C. BESSE, R. CARLES, N. MAUSER, AND H. P. STIMMING 



0.05E 




Figure 16. Blow-up time with varying a in quadratic os- 
cillations, supercritical power. Left: blow-up times for v. 
Right: comparison of v and w, dots with line for v, asterisks 
for w (Test 12). 



0.08 



0.07 - 




Figure 17. Blow-up time of v with varying a in quadratic 
oscillations, supercritical power (Test 13). 

energy level where the blow-up changes from two point blow-up to one point 
blow-up, and non-monotonicity can be observed. 



5. Numerical test, damped NLS 

5.1. Critical power. We now turn to ()2.1U|) . 

Test in one space dimension 

Test 16. For (|2.10j) in space dimension one, we first use the single 
Gaussian data (|3.H) . The scale C = 2 was chosen so that ||uo||^2 = 5.013, 
the smallest discretization parameters used are Np = 2^^ mesh points and 
At = 2.5 • 10~®. The blow-up time with respect to changing 6 is shown in 
Figure EDI The blow-up time is decreasing monotonically with 6, in accor- 
dance with the arguments of For 5 > 1.75, blow-up is prevented. 

Test 17. Next we use the "two-hump" data (|3.2j) . The data scale was 
chosen as C = 5 so that ||tio||^2 = 11-969. 

The finest discretization parameters used in this test are Np = 2^'^ mesh- 
points and At = 2.5 • lO"'^. 
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Figure 18. Blow-up time of v with varying a in quadratic 
oscillations, supercritical power, data with additional phase 
(Test 14). 




Figure 19. Blow-up time for w with varying a in quadratic 
oscillations, supercritical power, asymmetric data (asterisks 
and circles). Solid line: Blow-up time for v. (Test 15). 



The blow-up time with respect to changing S is shown in Figure 1211 It 
can be seen that the blow-up time is not monotonously increasing with 6. 
The effect is somehow more pronounced than in the case of Equation (|1.H) . 

Test in two space dimensions 

Test 18. For the two-dimensional case of (|2.1fl|) . we use the data (|8.8j) 
with C = 11, so that the initial mass is ||iio||^2 = 37.04. This is above 
the minimal value necessary for blow-up (^3). The finest discretization 
parameters used in this test are Np = 2^^ meshpoints and At = 5.0 • 10^®. 
The discretization domain is [—8,8]^. 

The blow-up time in dependence of 6 is shown in Figure TTH For 6 > 1.4, 
blow-up is prevented. 
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Figure 20. Blow-up time with varying damping constant 6 (Test 16). 



critical NLS with damping 



i 0.22 




0.1 0.2 0.3 0.4 0,5 0.6 0.7 0.8 0,9 1 




Figure 22. Blow-up time with varying damping constant 6 (Test 18). 
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supercritical NLS with damping 
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Figure 23. Blow-up time with varying damping constant 
6, supercritical power (Test 19). 

5.2. Supercritical Power. 

Test 19. We tested equation ([TTn|l in space dimension one with a = 
3. The data are chosen as (|3.2|) with C = 3.8, so ||'Uo 11/^2 = 3.53. The 
discretization parameters used in this test are Np = 2^^ meshpoints and 
At = 5.0 • 10"*^ or At = 10"^. The discretization domain is [-8,8]. 

The blow-up times with respect to the damping constant 5 are shown in 
Figure EHl Also in this case monotonicity is not true. For 5 > 1.17, blow- 
up is prevented. In this case, the occurrence of non-monotonicity is very 
sensitive to the size of the data, as already observed in Test 16. If we choose 
C = 3.6 or 4.0, instead of C = 3.8, we observe monotonicity. 

Test 20. We repeat Test 19 with a different scale of the data, we take 
C = 3.6 in (|3.2|) . so that ||mo||^2 = 3.17. The discretization parameters used 
in this test are Np = 2^^ mesh points and At = 5.0-10"^. The discretization 
domain is [—8, 8]. The blow-up times with respect to the damping constant 
6 are shown in Figure OH We can see that the blow-up is monotonous in S. 
The blow-up always happens a single point at the origin, and for 6 > 0.83, 
blow-up is prevented. 

All tests were done with both the TSSP and the Relaxation scheme (RS). 
The results agree, as an example we showed the comparison of the schemes 
in figure IT^ By using two numerical schemes with different discretization 
approaches, the possibility of observing just numerical defects introduced 
by a particular discretization method can be excluded. The effects of non- 
monotonicity found here can be seen with two inherently different numerical 
methods, which is an indication that the observations are not related to 
numerical errors but indeed analytical properties. 



In this paper, we have addressed numerically the question of the depen- 
dence of the blow-up time for solutions to nonlinear Schrodinger equations 
upon one specific parameter, in three cases: 



6. Conclusion 
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Figure 24. Blow-up time with varying damping constant 
6, supercritical power (Test 20). 

• Dependence upon the coupling constant A, for fixed n,a and uq: 

idtu + Au = -Xluf^u , {t,x) €R+ ; ^1^^0 = 1*0. 

• Dependence upon the magnitude of a quadratic oscillation intro- 
duced in the initial data: only a varies in the following equation, 

2 _ -Mi 

idtw + Aw = —\w\ '^w ; w\t=Q = UQ{x)e * . 

• Dependence upon the strength of the damping 5 > in 

idtilj + All) = - iSilj ; '4}\t=o = . 

In the LP' super-critical case cr > the tests we performed highly suggest 
that in either of the above three cases, the blow-up time is not monotonous 
with respect to the variation of the parameter considered. 

In the critical case cr = -, there is apparently no monotonicity in the 
first and in the third problem. In the second one however, our tests agree 
with the analytical result: there is monotonicity of the blow-up time with 
respect to a, as recalled in Proposition 12.41 

We used two numerical methods (time-splitting spectral method and a 
relaxation method), for the results observed to be more convincing. We 
may say that they are, since the two methods yield the same results (and 
not only just similar results, see Figure IT^ . 

Note that in some cases where we observed monotonicity reversal, the 
slope of the blow-up time/varying parameter curve may be rather steep near 
the monotonicity breakup. Compare Figure El near A = 2.7 with Figure ^ 
In the quadratic oscillations case. Figure shows a similar feature near 
a = 0.7. For the equation with a damping term, compare Figures OTI and 0^ 
with Figure I2UI This suggests that the underlying mechanism causing the 
monotonicity reversal is strongly nonlinear. It seems inappropriate to speak 
of instability though (for instance in Figure[Sl we saw that there is one point 
on the steep line after the local maximum). In these examples, the blow-up 
time changes rather fast compared to the modification of the L^-norm of 
the initial data (see also Remark |4.1() . Yet, the dependence of the blow-up 
time might also involve other analytical quantities, or geometrical features. 
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All the numerical counter-examples to monotonicity that we found contain 
a somehow nontrivial profile, inasmuch as the initial datum is formed of two 
(or more) humps. The question of an analytical justification remains open 
and challenging in these cases. 

Appendix A. Proofs 

A.l. Proof of Proposition 12.11 It is based on a fixed point argument 
using Strichartz estimates, which we first recall. 

Definition A.l. A pair {q, r) is admissible if 2 < r < (resp. 2 < r < 
oo if n = I, 2 < r < cxD if n = 2) and 

- = o r) := n 

Lemma A.l (Strichartz estimates, see e.g. |H2E1E])' Let {q,r), {qi,ri] 
and {q2,r2) be admissible pairs. Denote U{t) := e**^. 

1. There exists Cr such that for any (p G L^(M").- 

(A.l) \\U{-)'^\\L^{W-Lr) < CrlMl^ ■ 

2. There exists Cn^r-^ such that for any interval I and any 



(A.2) 



U{t - s)F{s)ds 

in{s<t} 



The proof of ProDOsition l2.1l relies on the following straightforward lemma: 

Lemma A.2. Let r = s = 2a + 2, and q = so that the pair {q,r) is 

admissible. Define k by 

^ _ 2a{2a + 2) 
~ 2-{n-2)a' 
Then k is finite, and the following algebraic identities hold: 
1 _ 1 2a 1 _ 1 2a 

r' r s ^ q' q k 

Proof of PropositionWl\ For A G {Id,Va.}, denote Ra := ||Aiio||i2. With 
the notations of Lemmas lA. II and IA.2I define: 

Xt := {u G C{%T[-H^) ; An G L°°([0, r[; L^) n L'?([0, T[; L''), 

P^IIl-{[0,T[;L2) < 2Ra , Pn||i,([o,T[;L'-) < 2^2, i^A, G {H, Vj} . 

Duhamel's formula for writes: 

(A.3) u{t) = U{t)uo + iX [ U{t- s){\u\^''u) {s)ds. 

Jo 

Denote F{u) the right hand side of ()A.3|) . We prove that F maps Xt into 
itself and is a contraction (in a weaker metric) provided that ()2.1() is satisfied, 
for e > sufficiently small. Denote := -L"([0, T[). From Lemmas I A. II 
and I A. 21 we have, for A G {Id,Va;}: 

P(^("))IIl-(l2) < \\Auoh2 + CX\\A{\u\^^u)\\ . , 
(A. 4) ^T\^ > 

< \\Auq\\^2 + C'A ||m||^1(^,) \\M\l'},{l-) ' 
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for some constant C depending only on n and a. From Gagliardo-Nirenberg 
inequality, we have: 

where 6{s) is given in Definition lA.ll We infer, for u £ Xt- 

\\A {F{u)) ||^.(^.) <R^ + CXT^'^/^Rlf-'^'^^Rlf'^RA , 

where from now on, we denote by C all the constants which do not depend 
on relevant parameters. Using Lemma I A. 11 we have similarly: 

for some other constant C depending only on n and a. We therefore have 
stability for 

(A.5) ^j.2a/kjf^^{l-Sis))j^2^a5is) ^ ^ 

We have contraction in the weaker metric L'^{[0, T[; U') under a similar con- 
dition: 

\\F{u2) - (Lr) < Cr,A II {\U2?"U2 - |ni ^'''U.i) 1 1 ^ 

< C\ + h2|lil(i.)) \\U2 - UiWL't^^Lr) 

<L\1 R^^ R^^ IF2 - (L-)- 

If (|A.5|) is satisfied, then F has a unique fixed point in Xt- Since (|A.5|) is 
nothing else but this proves the first part of Proposition 12.11 

When Mo S S, one can prove that the S-regularity is conserved along 
the time evolution by now considering A G {Id, V^.^, J(t)}, where J{t) = 
x + 2itVx is the Galilean operator. (This operator commutes with the group 
U{t) and acts on the nonlinearity we consider like a derivative.) Finally, we 
refer to |Hj for the conservation laws. □ 

A. 2. Proof of Corollary 12. In the above proof of local existence, we did 
not use the conservation laws. We now take them into account. 

Consider Duhamel's formula (|A.3|) . and seek an L^(L^ )-bound for V^^n. 
For A G {Id, Vx}, we have like above: 

Recall that s = 2a + 2; from the conservation of energy, 

< CX-'/' (||V,.n(t)||^/^ + \E\) < C {x-"'\\V Mt)\\%' + 1 
for A > 1, where C depends on uq and a, but not on A. We deduce, for any 

t e [o,r]: 



(A.6) 



Now assume that (|2.4() is satisfied for some e to be fixed later. Then by 
continuity (see Proposition 12. 1() . there exists r > such that 

(A.7) T^Hn-2)a fx^-\\V^u\\f^^L2) + l)< [2^" + 1) e. 



BLOW-UP TIME FOR NLS 25 

for < t < T. Choosing < e < £o{n,a), the last term of (|A.6|) can be 
absorbed by the left hand side for < t < r: 

Strichartz inequalities and the above estimate yield, for t £ [0, T] such that 
(TO) holds: 

By the conservation of mass, this estimate is interesting only when A = Vx- 
Up to choosing e even smaller, we see that as long as (|A.7|) is satisfied: 

||Va:n||^oo(-^2) < 2 ||Vi.tto||j;^2 . 

Now suppose that ()A.7|) is not satisfied for t = T: there exists a minimal 
time t* <T such that 

||Va;tt||^co(^2) < 2 ||V2;no||2,2 . 
The latest inequality yields, along with (|2.4j) : 

^2-(n-2)a J^^2a ||y^^||4<x^ ^^^^ ^ < ^2-(n-2)a ^X^- 2"'' \\V xUo\\t2 + l) < 24'^£. 

This contradicts the definition of t* . Therefore, ()A.7|) is satisfied for t = T, 
which completes the proof of Corollarv 12.21 

A. 3. Proof of Proposition 12. 4L For a G M, define: 

Ja{t) =x + 2i{t - a)Vx = 2i{t - a)e^(^Vx i e~'^(^ 

Denote 

/ J. \ na—2 

m = h{tr-^ = 1 - - 



Besides the conservations of mass for u and f , and of energy for n, we have 
the following evolution laws: 

^ fuv7 „,/+M|2 4>{t) |L,^^M|2<7+2 ^ _ -'^'(*)|L,^^M|2<7+2 . 



■(iiv..(t)iii2 - ^ib(t)iii^^i) = ^um] 



'-{t - a)mHmL2t^. - -—rit - afcp'mv{t)\\l^+^. 



Lemma A. 3. Let uq G S, 2/n < a < 2/{n - 2) and T > 0. The following 
assertions are equivalent: 

(1) u blows up at time T. 

(2) \\Vxu{t)\\L2 +00 ast^T. 

(3) ||u(t)||j^2<T+2 +00 as t ^ T. 
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+00 as t ^ T. 



(4) For any a^T, \\Ja{t)u\\]^2 

Proof of Lemma \A . S[ The equivalence between the first two points foUows 
from Proposition 12.11 The second point is equivalent to the third one by 
conservation of the energy. 

For a ^ T, Gagliardo-Nirenberg inequalities and the second writing for 
Ja yield: 

C{n,a) 



\um 



2a+2 

L2a + 2 



< 



t - a " 



2-(n-2)(T I 



Ja{t)u\ 



L2 ■ 



We infer (3)^(4) by the conservation of mass. Finally, if ||u(t)||^2cT+2 remains 
bounded as t — > T, then integrating the first evolution law above shows that 
II jQ,(t)ii 1 1 2,2 remains bounded as t ^ T. □ 

A similar result holds for v. Notice that if a > 0, then for < t < a, 
(f) is smooth (bounded) and does not cancel. If a < 0, the same holds in a 
compact of M+ . 

Lemma A. 4. Let uq G S, 2/n < a < 2/{n - 2) and T > 0, with T < a if 
a > 0. The following assertions are equivalent: 

(1) V blows up at time T. 

(2) ||Va;u(t)||L2 +00 ast^T. 

(3) ||u(t) 112,2.7+2 +00 as t ^ T. 

(4) For any a^T, \\Ja{t)v\\i2 +00 as t ^ T . 

Proof. The equivalence between (1) and (2) follows like for Proposition 
since (j) is bounded from above and from below by positive constants. 

The other equivalences follow like above, thanks to Gagliardo-Nirenberg 
inequalities. □ 

Then Proposition 12.41 follows from Lemma lA.41 and the identity: 

at 



\Ja{t)v\\L'2 



t 



L2 



A. 4. Proof of Proposition 12. 6L It consists of a slight adaptation of the 
proof of Proposition 12. II Considering (|2.1()|) . the analog of ()A.4|) yields: 



\A{F{u) 



< \\Au, 

< \\Au, 



oIIl2 



oIIl2 



+ CA 



+ CA 



-2crSt 



2(7 



L« (L-') 



where the notation now stands for LP{ 
inequality, we have: 



P^IIl«(L'-) ' 

Using Gagliardo-Nirenberg 





< C 











u 



,1-6(3) 



IV.xU 



Ms) 



\L^{L^) ■ 



Mimicking the proof of ProDosition l2.11 we can apply a fixed point argument 
provided that: 



-St 



l-5{s)j^5{s) 



<ci. 



This property is a consequence of the stronger one: 



R 



Id 
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which is the condition stated in Proposition 12.61 whose proof is now com- 
plete. 

Remark. This proof yields similar results for a larger class of damping. For 
6 > and F E L^(M_|_) a positive function, define: 



H{t,5) = -6- 



Then for 6 satisfying (|2.11j) with another constant C (depending of F), the 
solution of 

does not blow up in the future. In Prop. ESI we considered the case F{y) = 
e~y . Considering F{y) = (1 + y)~" with a > 1 yields a damping function 

which goes to zero as t ^ +00. Even though the damping is less and less 
strong, (|2.11jl ensures global existence in in the future. 
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